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ABSTRACT 

We  give  a  proof  that  Total  Variation  (Uq(«)}  <<  1  can  be  replaced  by 
Sup{uQ( • ) }  <<  1  in  Glimm's  method  whenever  a  coordinate  system  of  Riemann 
invariants  is  present.  The  argument  is  somewhat  simpler  but  in  the  same 
spirit  as  that  given  by  Glimm  in  his  celebrated  paper  of  1965. 
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SIGNIFICANCE  AND  EXPLANATION 

i 

A  system  of  two  conservation  laws  in  one  dimension  is  a  set  of  first 

Ctrl, i  ^ 

order  nonlinear  partial  differential  equations  of  -the  form- 

_  _ .  —  -  '  u*  +  f(u,v)„  -  0  , 

(Dl  ' 

vt  +  g(u,v)  -  0  , 

tfs. 

where  (u,v)  is  a  vector  function  of  (x,t),  x  ^  R,  t  f  0.  The  Cauchy 
problem  asks  for  a  solution  of  (1)  given  the  "initial*  values  of  u  and  v 
at  time  t  *  0.  Equations  of  type  (1)  arise,  for  example,  in  gas  dynamics 
where  they  express  the  conservation  of  quantities  like  mass,  momentum  and 
energy,  when  diffusion  is  neglected.  Typically,  smooth  solutions  of  (1) 
cannot  be  found.  This  is  due  to  the  formation  of  shock  waves.  Shock  waves 
are  the  mechanism  by  which  entropy  is  dissipated  in  solutions  of  ( 1 ) .  This 
paper  gives  a  proof  that  solutions  exist  even  after  shock  waves  form,  so  long 


as  the  amplitude  of  the  waves  are  not  too  great  initially. 
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Blake  Temple 

We  consider  the  Cauchy  problem 

ut  +  F(u)x  -  0  (1) 

u(x,0)  -  Uq( x)  ,  (2) 

where  (1)  denotes  a  strictly  hyperbolic  system  of  two  conservation  laws,  u  -  (uj,u2), 

F  -  (f,g).  Let  ^p'Rp'  P  “  be  the  eigenvalues  and  corresponding  eigenvector  fields 

associated  with  the  matrix  VF,  X1  <  X^.  Assume  that  U  is  a  neighborhood  of  a  state  u 
in  u-space  in  which  each  characteristic  field  is  either  genuinely  nonlinear  ( V >  0) 
or  else  linearly  degenerate  (7X^»R^  50),  and  such  that  X^(u)  <  X2(v)  for  all 
u,v  t  U  [4].  Without  loss  of  generality,  assume  u  -  0.  In  this  note  we  give  a  simplified 
proof  of  the  following  result  which  is  contained  in  the  results  of  Glimm  [2]  and  which  is 
required  for  the  proof  in  [12].  (A  stronger  result  also  follows  from  the  analysis  in  [3] 
which,  however,  involves  the  theory  of  approximate  characteristics  and  is  much  more 
technical.)  Let  u(x,t)  denote  a  weak  solution  of  (1),  (2)  which  is  a  limit  of 
approximate  solutions  generated  by  the  random  choice  method  of  Glimm. 

Theorem  Is  For  every  VQ  >  0  there  exists  a  small  constant  6  >  0  and  a  large  constant 
G  >  0  such  that,  if 

TV{u0(*)>  <  V0  ,  (3) 

and 

lu0(.)ls  <  «  ,  (4) 

then 


for  all  t  >  0.  Here 


lu<»,t)l$  <  Glu0(.)lg  , 
Iu0<*)ls  =  Sup  {Uq(  • ) }  . 


(5) 

(6) 
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It  suffices  to  verify  Theorem  1  for  any  approximate  solution  uh  generated  by  the 
random  choice  method.  Recall  that  there  exists  a  coordinate  system  of  Riemann  invariants 
for  (1)  in  a  neighborhood  of  u  •  0.  (Indeed,  Theorem  1  and  the  proof  to  follow  apply  to 
any  system  satisfying  the  above  assumptions  and  for  which  there  exists  a  coordinate  system 
of  Riemann  invariants  [2].)  Moreover,  the  Riemann  problem  is  uniquely  solvable  in  a 
neighborhood  of  u  ■  0  by  the  method  of  Lax  [4].  Such  a  solution  consists  of  a  1-wave 
Y 1  followed  by  a  2-wave  Y^  each  of  which  is  either  a  shock  wave  or  a  rarefaction 
wave.  Assume  that  U  is  a  neighborhood  of  u  “  0  satisfying  all  of  the  above 
conditions.  Define  the  strength  of  a  wave  |yPI  to  be  the  absolute  value  of  the  change  in 
the  opposite  Riemann  invariant  between  the  left  and  right  states  of  the  wave.  Finally,  in 
order  to  set  notation,  we  briefly  review  the  construction  of  the  random  choice  method 
approximate  uh . 

Let  h  be  a  mesh  length  in  x,  and  let 

k  -  Ch 

be  the  corresponding  mesh  length  in  t,  C  >  Sup  { | X  < u) | } .  For  i,j  t  t,  J  >  0,  let 

p,u«U  P 

5  ih,  tj  =  jk.  Let  a  be  a  sample  sequence,  a  =  {a^}^,  0  <  Sj  <  1.  For  given 

initial  data  uQ( • )  C  U,  define  the  random  choice  method  approximate  solution 
h  h 

un(x,t)  =  u  (x,t;a)  by  induction  on  j  as  follows:  First,  for  xi  <  x  <  xi+1,  define 

uh(x,0)  =  Ug(x)  -  u0(xi  +  |)  . 

Next,  assume  for  induction  that  uh(x,t)  has  been  defined  for  t  <  tj.  Define 

uNx.tj)  =  uh(Xi  +  *jh,t j-)  , 

and  for  tj  <  t  <  tj+1,  define  uh(x,t)  to  be  the  solution  of  the  Riemann  problem  posed 
in  (3.3)  at  time  tj.  By  (3.1),  uh  is  well  defined  so  long  as  uh(x,tj)  C  U  for  all 
tj.  Let  YPj  denote  the  p-wave  that  appears  in  the  solution  of  the  Riemann  problem  that 
is  posed  at  (x^,tj)  in  the  approximate  solution  uh.  Recall  that  the  quadratic 
functional  associated  with  u^  is  defined  by 

Q(t)  i  l  Iy^IIy^I  (7) 

where  the  sum  is  over  all  waves  that  approach  at  time  tj,  tj  <  t  <  tj+1«  Let  A^j  denote 
the  interaction  diamond  centered  at  (x^ftj),  and  let  D^j 


denote  the  products  of 


A^[12].  Vie  use  the  following  notation: 

V.  s 

3 

'^1  ' 

P#i 

(8) 

Qj  5 

Qttj)  »  Qttj+J  , 

(9) 

Dj  = 

I  Dij 

i 

(10) 

Fj  5 

vj +  ®j 

(11) 

Sj  = 

luh(*,t.)l.  . 

3  ® 

(12) 

Note  that  V  estimates  the  total  variation  of  u  ( • » t ^ )  and  that  (7),  (9)  give 


immediately  that 


(13) 


We  show  that  Theorem  1  is  a  consequence  of  the  following  lemma  which  is  a  restatement  of 
results  in  [2]: 

Lemma  1 :  There  exists  a  constant  Gg  >  1  depending  only  on  F  such  that,  if 
uh(x,t)  e  U  for  all  t  <  tj,  then 

Vi  -  V)  <  G0SjDj  ' 


Vi  ~  Sj  ‘  G0SjDj  ' 
Qj+i  -  Oj  <  <Gosjv: j  -  1}Dj 


(14) 

(15) 

(16) 


Proof  of  Theorem  1:  Fix  Vn  >  1.  Choose 


.  2  “t 
Sq  <  (“Vn*  )  *  5  ' 


where  Gn  is  large  enough  so  that 


{u  i  |u|  <  Gg1)  C  U 
We  show  by  induction  that  (17)  implies 


(17) 


(18) 


Sj  <  e 


2G 

0  0_ 


s0  ' 


and 


OQSjVj  <  4 


(19) 


(20) 
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k/V 


Also  not*  that  (10)  and  (11) 


for  all  j  >  0.  Not*  that  (19)  givaa  (5)  with  6  -  a 
loply 


2G.V? 
0  0 


G0Sj  <  j  • 


1 


Thua,  ainca  (8)  gives 


(21) 

(22) 

(22) 

(24) 

(25) 


Cj-1  -  «  <°0SjVj  "  ^Dj  ' 

pJ+i  -  pj  «  {0o8j  +  Vjvj  -  1}DJ  ' 
estimates  (19)  and  (20)  alao  imply  that 

fij+1  *  Qj  «  ”  i  °j  ' 

Pj+1  -  pj  *  "  7  DJ  • 

He  now  verify  (19)  and  (20)  by  induction.  The  idaa  her*  ia  that  (20)  guarantaaa  that  both 
(Qj)  and  {r_j }  era  decreasing.  The  decreasing  of  Qj  controls  Sj  at  tha  Induction 
step,  while  the  decreasing  of  {P^}  maintains  (20)  at  the  induction  step,  sine*  than 
Vj  <  <  V 

First,  when  j  -  0,  (10)  and  (11)  follow  from  (17).  So  asstaae  (10),  (11)  hold  for 

j'  <  j.  He  verify  (10),  (11)  for  j*  -  J  ♦  1.  8y  (24)  and  (15), 

^♦1  "  ®k  <  -  W  ' 


or 


®k-M  <  0  ♦  »olOk- W>*k  ' 


(26) 


for  k  <  j  +  1.  Thus  by  (26) 


8j+1  <  n  {1  ♦  2C0(Qk  -  Qlt^.1J>s0  . 
k“0 


(27) 


But  one  can  easily  verify  that  the  mxImuh  of 


n  (1  ♦  a*) 

k-0 


over  all  nonnegative  sequences  {a^}-^  Mti**yi"9  f  a^  <  M  ia  attained  whan 

k»0 


*»  “  jTT 


for  all  k.  Thus  by  (13)  and  (24), 
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'll 

yj 

yJ 

,:.(i 

*i 


M 

■*5 


l  ICk  -  W  «  O0  <  vo  ' 

k»0 


j  2G  V2  ^+1  28  V2 

kl0  {1  +  4  i1  +r?-T}  4  e  0  0 


Therefore  by  (27), 


sj+1  *  e  S0  ' 


which  verifies  (19)  at  )  *  I.  Moreover,  (25)  implies 


Vj+i  <  Fj+1  <  F0  <  2V0  . 


s:l 


Thus  by  ( 17) , 


2G  V 

0  0  2  1 

G0Sj-MVj+1  <  2Q"  V0  <  4  ' 


which  verifies  (20)  at  j  +  1.  This  completes  the  proof  of  Theorem  1. 


’  v '■ , *  >  y9**?*.  ~ 


REFERENCES 


[1]  Courant,  C.  and  Friedrichs,  K.  O.,  "Supersonic  Flow  and  Shock  Haves", 

Springer-Verlag,  1948. 

[2]  Glimm,  J. ,  Solutions  in  the  large  for  nonlinear  hyperbolic  systems  of  equations. 

Comm.  Pure  Appl.  Math.,  18  (1965),  695-715. 

[3]  Glimm,  J.  and  Lax,  P. ,  Decay  of  solutions  of  systems  of  nonlinear  hyperbolic 
conservation  laws,  Amer.  Math.  Soc.  Memoirs,  101  (1970). 

[4]  Lax,  P.,  Hyperbolic  systems  of  conservation  laws,  II,  Conan.  Pure  Appl.  Math.,  10 
(1957),  537-566. 

[5]  Luskin,  L.  and  Temple,  J.  B. ,  The  existence  of  a  global  weak  solution  to  the 
nonlinear  waterhammer  problem.  Comm.  Pure  Appl.  Math.,  35  (1982),  697-735. 

[6]  Smoller,  J.  A.,  Shock  Haves  and  Reaction-Diffusion  Equations,  Springer-Verlag,  1980. 

[7]  Temple,  J.  B.,  Global  solution  of  the  Cauchy  problem  for  a  class  of  2  x  2  nonstrictly 
hyperbolic  conservation  laws,  Adv.  Appl.  Math.,  3  (1982),  335-375. 

[8]  _ ,  Solutions  in  the  large  for  the  nonlinear  hyperbolic  conservation  laws 

of  gas  dynamics,  J.  Diff.  Eqns.,  Vol.  41,  No.  1  (1981),  96-161. 

[9]  _ ,  Systems  of  conservation  laws  with  coinciding  shock  and  rarefaction 

curves.  Contemporary  Mathematics,  Vol.  17  (1983),  143-151. 

[10]  _ ,  Systems  of  conservation  laws  with  invariant  submanifolds,  Trans.  Am. 

Math.  Soc.,  Vol.  280,  No.  2  (1983),  781-795. 

[11]  _ ,  No  L^contractive  metrics  for  systems  of  conservation  laws,  Trans.  Am. 

Math.  Soc.,  Vol.  288,  No.  2,  April  1985. 

[12]  _ ,  Decay  with  a  rate  for  non-compautly  supported  solutions  of 

conservation  laws,  MRC  preprint. 


security  classification  of  this  page  rmiMi  dx«  &■«•»« 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 

#2855 


4.  TITLE  (md  Submit) 

SUP-NORM  ESTIMATES  IN  GLIMM'S  METHOD 


7.  AUTHORfc; 

Blake  Temple 


PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

(See  Item  18  below) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J.  RECIPIENT’S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERS) 

DAAG29-80-C-00  41 
DMS-8210950,  Mod.  1. 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

Work  Unit  Number  1  - 
Applied  Analysis 


12.  REPORT  DATE 

August  1985 


MONITORING  AGENCY  NAME  •  ADDRESSfi/  dllloront  from  Controlling  OlHeo)  15.  SECURITY  CLASS,  (of  chi*  roport) 

UNCLASSIFIED 

I  IS*.  DECLASSIFICATION/DOWNGRADING 
I  SCHEDULE 


18.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office 

National  Science  Foundation 

P.  0.  Box  12211 

Washington,  DC  20550 

Research  Triangle  Park 

North  Carolina  27709 

•9.  KEY  WORDS  (Continue  on  reretme  eld •  It  necoooory  and  Identity  b y  block  number) 

Riemann  problem,  Randan  choice  method.  Stability,  Conservation  Laws 
Cauchy  problem 


20.  ABSTRACT  ( Continue  on  rovoroo  old*  II  noeotoory  and  Idonllty  by  block  number; 

We  give  a  proof  that  Total  Variation  {uQ(»)}  <<  1  can  be  replaced  by 
Sup{uQ(*)}  <<  1  in  Gibran's  method  whenever  a  coordinate  system  of  Riemann 
invariants  is  present.  The  argument  is  somewhat  simpler  but  in  the  same 
spirit  as  that  given  by  Glimm  in  his  celebrated  paper  of  1965. 


